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1 Statement of the problem 



Let M be a smooth 2-manifold and let C a smooth cross field defined on M. For all general definitions about 
cross fields, we refer to G). 

We assume that C has only a finite number of isolated singularities, and that all such singularities are of the first 
order, i.e., the index of a singularity can be either |, or — |. Each point of M that is not a singularity of C will be 
said to be regular. 

A streamline (a.k.a. integral line) of C is a smooth line on M that is tangent/orthogonal to the directions defined 
by C at each point; a streamline either is a cycle, or it has its endpoints at (not necessarily distinct) singularities. A 
streamline that ends at singularities is called a separatrix. 

Let S be the set of singularities of M. Under the assumptions above, there exist a finite number of separatrices 
that join pairs of points of S, which may cross orthogonally at a finite set R of regular points of M. Separatrices 
partition M into a finite number of "quadrangular" domains, where by quadrangular, we mean that each region is 
bounded by four segments of separatrix. 

Let E be the set of segments of separatrix that join two points of V = S U R. Then G = ( V, E) is a graph 
embedded on M with the following properties: 

1. Each minimal cycle is formed by four edges; 

2. Each regular vertex in R has valence four; 

3. Each singularity in S having index \ has valence three; 

4. Each singularity in S having index — \ has valence five. 

Note that, since G is embedded on M, the arcs incident at any vertex v E V can be arranged in radial order 
around v. Let c = (vq, . . . , v n ) be a chain in G, and let V{ be a regular vertex on the chain. We say that c crosses 
Vi if the edges (i^-i, vi) and (yi, Vi+i) are not consecutive in the radial order around otherwise we say that c 
turns at V{. A separatrix corresponds to a chain (so, n,..., r^-i, Sk) where so and Sk are singularities, all the r\ 
are regular and the chain crosses at all regular vertices. 

We aim at obtaining another graph G' = (V f , E') such that V' = S U R f , with R' possibly empty, which fulfills 
all properties 1-4 defined above (with R r in the place of R) and such that the lines of E' are almost aligned with C. 
The problem can be defined as an optimization one, by giving an energy £{G') that is proportional to the number 
of vertices of R' and to the "warping" introduced by G' on C, which depends on how much the edges of E' deviate 
from lying on streamlines of the cross field. 

In the following section, we remain generic about energy £, as well as about geometry of lines of E, while we 
describe a greedy algorithm that provides a heuristic for tackling with this optimization problem in general, by 
treating it as a graph simplification problem. 
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Figure 1. Sub-atomic operation used to repair a S-defect. Circles and squares represent singular 
and regular vertices; defected vertices are depicted in orange and yellow; defects are represented as 
red segments; the defect at vertex v can be repaired by cutting one of the two separatrices depicted 
in black. 



2 Graph simplification 

We get in input graph G as defined in the previous section, and we generate a graph G' = (5 U i?', that 
fulfills properties 1-4, trying to reduce the number of regular vertices in R f , while disregarding the geometry of 
edges in E' . 

Our algorithm works through macro-operations that allow us editing the graph while maintaining it consistent 
with properties 1-4. Each macro-operation corresponds to a series of atomic operations that may either delete, or 
detach and reattach edges in the graph, and possibly remove regular vertices. 

Let us consider a vertex v together with its incident edges. If we detach an edge from v, we generate a defect 
at v. If v is either singular, or regular, then we call this defect either a S-defect, or a R-defect, respectively. If we 
attach an edge e to a defected vertex v, then we repair a defect. Note that, since edges are arranged radially, each 
defect has a given position in the radial order about v. Repair is valid only if the position of edge e in the radial 
order, which depends on the embedding of e on M, is coincident with the position of the repaired defect. Note 
that we are interested in such an embedding only just from the point of view of topology, i.e., the shape of e is not 
relevant, as long as it does not cross other edges. 

Repairing a defect corresponds to the sub-atomic operation described in the following. In our algorithm, there 
can be just one defect per vertex (actually, there will be just one S-defect and one R-defect in the whole graph 
after each atomic operation). A defect at vertex v is repaired by cutting one of the two separatrices that run 
"parallel" to the defected edge, and attaching one of the two branches of the cut separatrix to v at the defect. This 
is accomplished as depicted in Figure[TJ let e f the edge adjacent to the defect, in the radial counterclockwise order 
about v, and let v' the other endpoint of e'\ if v' is regular, let e" the edge next to e' in counterclockwise order 
about v f , and let v" be the other endpoint of e"\ edge e" is deleted and v is connected to v" with a new edge; 
this generates a R-defect at v f , and the branch of cut separatrix ending at such defect will be called the dangling 
branch. 

The same operation could be done symmetrically by rotating clockwise about v and cutting the other "parallel" 
separatrix, therefore there are two possible options for repairing a defect. If a neighboring vertex v' is singular, the 
corresponding option cannot be used to repair the defect; if both neighboring vertices are singular, then the defect 
cannot be repaired, and the algorithm will have to backtrack (see next about backtracking). 
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Figure 2. The four possible configurations generated by the Delete-separatrix operation. 

Atomic operations are of the following two types: 

1. Delete-separatrix (Figure [2]): delete all edges and regular vertices of a separatrix; for each regular vertex v 
deleted, the two other edges incident at v are merged to form a single edge; one of the two S-defects caused 
by removing the separatrix is repaired (as explained above), by generating a R-defect; 

2. Switch separarix (Figure [3}: repair a R-defect by cutting a separatrix (as explained above) and delete the 
dangling branch of the separatrix; delete all regular vertices traversed by the dangling branch and merge 
pairs of edges incident at them, as in the previous operation; if the other end of the dangling branch is a 
singular vertex, then a S -defect is generated. 

A macro-operation consists of one Delete-separatrix operation followed by a loop of Switch-separatrix op- 
erations, which ends when no S-defect remains in the graph. Note that after each operation, one S-defect and one 
R-defect are present in the whole graph. Each Switch-separatrix removes a S-defect and potentially produce an- 
other S-defect; however, if the dangling branch of separatrix removed by this operation ends at the (only) R-defect, 
then no defects are left. This sets the stop condition for the cycle in a macro-operation. 

Now let us consider the possible choices we have in applying a single macro-operation. Let n be the total 
number of separatrices of field C (note that this number remains constant throughout processing). We thus have n 
options to give a separatrix in input to Delete-separatrix, and four options to implement this operation, as shown 
in Figure |2| hence a total of 4n options to initialize the simplification process. Given one such choice, once we 
enter the loop, we have two options for each Switch-separatrix operation, as shown in Figure |3] This gives a total 
space of solutions consisting of a forest of 4n binary trees, one for each initial choice. Leaves of a tree correspond 
to either termination a solution, or a stuck condition, which occurs when a defect cannot be repaired. All trees 
are guaranteed to have a finite number of nodes, because each Switch-separatrix operation removes at least one 
regular vertex and there exist a finite number of such vertices. 

An algorithm searching for the best simplification (with respect to some given criterion) could explore this 
space of solutions and select the best one (i.e., the one that minimized the given energy E). The complexity is 
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Figure 3. The two possible configurations generated by a Switch-separatrix operation. 

combinatorial, but it could be feasible in practice, because separatrices are usually not very numerous, and the 
depth of trees is also usually quite small. 

Our approach aims at evaluating a suboptimal solution with a greedy algorithm that explores just a path in a 
tree from the forest, possibly backtracking and choosing an alternative path if it gets stuck at a leaf. The separatrix 
to be deleted and the way to repair the first S-defect are selected with some heuristic; at each cycle, the way to 
repair the R-defect is selected with some other heuristic. The algorithm keeps track of the operations performed 
in a stack. If it gets stuck at a leaf, it backtracks by removing the last operation from the stack, and by freezing the 
corresponding option, until a node with an open (non-frozen) child is found. In the worst case, the whole tree is 
visited while no path leads to a solution, so the graph is left unchanged. 

The repeated application of a macro-operation simplifies the graph progressively. This process can be repeated 
until some used-defined criterion occurs, such as having reached either a given number/percentage of regular 
vertices, or a maximum tolerated level of warping with respect to field C. 

3 Optimization criteria 

We now consider the geometric aspects involved in the simplification process, in order to devise some heuristics 
for driving our algorithm. We need a criterion to choose the first separatrix and S-defect to be repaired, and a 
criterion to chose between the two possible options to repair a given defect. 

Since our aim is to reduce the number of regular vertices of G, we start the algorithm by deleting the separatrix 
that crosses the largest number of vertices. The criterion for selecting the edge to detach to repair a defect is the 
same in both atomic operations and it is based on the concept of minimal drift explained in the following and 
depicted in Figure [4] 

Let us consider a defect at a vertex v, and let v f , e' and e" be described as in the previous section. Let so be 
the separatrix leaving v in the direction of the defect and let s f the separatrix leaving v in the direction of e', both 
separatrices being computed on the original graph G. Similarly, let s\ be the branch of separatrix through v' in 
the direction of e", let v s be the singularity at the end of such branch, and let s[ be the separatrix leaving v s in the 
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Figure 4. The minimal drift is the ration between the area of highlighted in blue and the square of 
the length of the dotted diagonal. The yellow square denotes regular vertexes in the graph that has 
been removed by simplification operators. 

direction orthogonal to si. The four separatrices so, <Sq, si and s[ form a quadrangular region Q on M. Let Aq be 
the area of Q, and let s be the shortest line through Q joining v to v s Finally, let l§ be the length of s. We define 
the drift of s as the ratio between the area of Q and the square length of s, 

Aq 

The drift is roughly proportional to the ratio between the "thickness" and the "length" of strip Q and it roughly 
measures how much line s deviates from being a streamline of C. 

Assuming that we can measure this drift efficiently, we choose each time to repair a defect by using the edge 
corresponding to the line with minimum drift. 

4 Implementation 

Our prototype implementation works on cross-fields induced by quadrangular meshes. Let M be a quadran- 
gular mesh, such that all its vertices have valencies between three and five. Mesh M stands here as a discrete 
approximation of a smooth manifold. We may assume a (discrete approximation of) a cross field C to be defined 
on M by assuming that: 

• The edges of M are aligned with the streamlines of C; 

• C is regular inside each face of M. 

Under these assumptions, the singularities of C correspond to the irregular vertices of M, and separatrices corre- 
spond to chains of edges that start at those irregular vertices and cross at all regular vertices. The input graph G 
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is obtained easily by traversing each separatrix and tagging all its edges and vertices. Thus, G will consist of all 
tagged vertices and edges of M after traversing all separatrices. 

In Figure [5} we show the graph of separatrices obtained from the quadrangulation of a popular dataset, obtained 
with the method presented in [lj. Each separatrix is depicted with a different color. Note how some separatrices 
may take involved routes. In particular, some of them form long helices, which generate many intersections with 
other separatrices, thus many regular vertices of G. These separatrices will be the first ones to be selected for 
deletion by our algorithm. 

If the input mesh M is highly regular, as in the case of meshes obtained with [lj, the algorithm can be greatly 
simplified by assuming that all edges of M have the same (unit) length. Note that each edge of G can consist of a 
chain of edges of M, thus it will have an integer length. Under these assumptions it is trivial to compute both the 
length of a separatrix and the drift, as described in the previous section. 

Once the graph has been simplified, we actually compute each new separatrix s as depicted in Figure |4j by 
assuming it as a "diagonal" of strip Q: in fact, strip Q is a regular rectangular sub-grid of M, with the endpoints 
of s at two opposite corners of the grid. 

5 Results 

The first experimental results obtained with our method are depicted in Figure [5] The input mesh Fertility con- 
sists of 3357 quads and it has been obtained with the remeshing method in (T). The mesh contains 48 singularities. 
The graph G of separatrices induced by such mesh contains 108 separatrices consisting of 2217 regular vertices. 
Among them, there are some long helices, which generate many intersections with other separatrices, thus many 
regular vertices of G. These separatrices will be the first ones to be selected for deletion by our algorithm. 

Figure [5] (a) show the graphs of separatrices of the original mesh. After one step of simplification the graph 
contains only 197 regular vertices and it shown in Figure [5] (b). A second simplification step allows to reduce the 
number of regular vertices to 88. It is not possible to further simplify the mesh without introducing a very high 
drift. 

Since our algorithms minimize the drift during the simplification, singularities are naturally aligned as it can be 
seen in the base or in the back of the mother. 

6 Concluding remarks 

This draft paper presents just preliminary work. Our method can work equally on line fields containing sin- 
gularities of the first order of type lemon and star. We plan to extend our method to work on line and cross 
fields computed on generic meshes. There exist many methods for computing such fields, most of which work on 
triangle meshes. A challenging problem is to extract graph G of such fields without computing a quadrangular 
remeshing aligned with the cross field. 
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Figure 5. The graph of separatrices computed on Fertility model remeshed with [1] (a). The graph 
simplified with a single step of simplification (b) . A second simplification step reduces the number 
of regular vertices with respect to the original model of about 96% (c). 
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